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ABSTRACT 

We investigate the imprint of nonlinear matter condensations on the Cosmic Mi- 
crowave Background (CMB) in an = 1, Cold Dark Matter (CDM) model universe. 
Temperature anisotropics are obtained by numerically evolving matter inhomogeneities 
and CMB photons from the beginning of decoupling until the present epoch. The under- 
lying density field produced by the inhomogeneities is followed from the linear, through 
the weakly clustered, into the fully nonlinear regime. We concentrate on CMB tem- 
perature distortions arising from variations in the gravitational potentials of nonlinear 
structures. We find two sources of temperature fiuctuations produced by time- varying 
potentials: (1) anisotropics due to intrinsic changes in the gravitational potentials of 
the inhomogeneities and (2) anisotropics generated by the peculiar, bulk motion of the 
structures across the microwave sky. Both effects generate CMB anisotropics in the 
range of fO"'' ^ AT/T ^ fO~^ on scales of ~ f°. For isolated structures, anisotropics 
due to proper motion exhibit a dipole-like signature in the CMB sky that in principle 
could yield information on the transverse velocity of the structures. 

Subject headings: cosmic microwave background - cosmology 



e-mail: rtuluie@astro.psu.edu 
e-mail: pablo@astro.psu.edu 



1 



Recently, there has been a renewed interest on 
studying anisotropies of the CMB produced by the 
late time evolution of density inhomogeneities 
(Martinez-Gonzalez and Sanz 1990; Martinez-Gonzalez 
et al. 1990; Martinez-Gonzalez et al. 1992; Panek 
1992; Fang and Wu 1993; Meszaros 1994; Martmez- 
Gonza.lez et al. 1994; Hu and Sugiyama 1994;). How- 
ever, most of these and other studies in the past (Rees 
and Sciama 1968; Dyer 1976; Kaiser 1982; Thomp- 
son and Vishniac 1987; Dyer and Ip 1988; Nottale 
1984) have only considered the effects on the CMB by 
a single, spherically symmetric structure using Swiss 
cheese models, thin-shell approximations or Tolman- 
Bondi solutions. The imprint on the CMB of mul- 
tiple, non-trivially distributed matter structures, as 
well as the unambiguous linkage between the linear, 
weakly clustered and fully nonlinear epochs, has not 
been yet taken into consideration; in particular, the 
peculiar motion of matter structures has been ignored 
so far. 

In this letter, we address the signatures in the 
microwave sky from time-varying gravitational po- 
tentials associated with nonlinear density inhomo- 
geneities in the universe, first studied by Rees and 
Sciama in 1968. This effect, hereafter referred as the 
Rees-Sciama effect, is a secondary source of temper- 
ature anisotropy; however, it could yield important 
information on the underlying density field that gen- 
erates this effect. Our goal is to study the contri- 
bution of the Rees-Sciama effect to the total CMB 
anisotropy, in particular, finding the correlation be- 
tween the Rees-Sciama effect with the underlying den- 
sity field and proper motion of structures that pho- 
tons encounter along their paths. For this reason, we 
have ignored Doppler shift anisotropies from Thom- 
son scattering. Doppler shift anisotropies from late 
Thomson scattering are important if a substantial 
reionization of the IGM occurred (Tuluie et al. 1995). 

Our study is based on numerical simulations of un- 
biased, SI = 1, /} = 0.75, CDM models of the uni- 
verse. For the parameters of our models, voids are 
typically of size ~ 30/i"i Mpc, and structures pro- 
ducing the stronger signatures ~ 10h~^ Mpc. Our 
study follows the approach developed by Anninos et 
al. 1991 that relates density inhomogeneities to CMB 
fiuctuations by evolving CMB photons through the 
forming matter inhomogeneities. It is important to 
emphasize that no assumptions are made about the 
nonlinear nature of the density perturbations or pecu- 
liar bulk motions. We only require sub-horizon scales 



and non-relativistic proper velocities, so linear theory 
of cosmological perturbations holds (Mukhanov et al. 
1992). Furthermore, since we propagate CMB pho- 
tons through matter that evolves from the beginning 
of decoupling up to present time, our results natu- 
rally take into consideration not only the gravitational 
lensing of the CMB caused by the late, nonlinear den- 
sity perturbations, but also the cumulative effects on 
the CMB photons as they propagate through multiple 
and different regions of voids and clusters. 

In linear theory of cosmological perturbations 
(Mukhanov et al. 1992), the evolution of the density 
field is governed by a gravitational field (f) which rep- 
resents small perturbations against the homogeneous 
background ((f) ^ 1). The temperature anisotropies 
of the CMB associated with the gravitational poten- 
tial (f) are given by (Sachs and Wolfe 1967) 

f = ^.-^. + 2f'Pt (1) 

where the integral is taken along photon trajecto- 
ries. The first and second terms in equation (1) repre- 
sent the Sachs- Wolfe and Rees-Sciama effect, respec- 
tively. The numerical factors multipling each term in 
equation (1) differ among various authors depending 
on the choice of gauge (Padmanabhan 1993). Our re- 
sults were obtained using the longitudinal gauge only. 
In an SI = 1 matter-dominated universe, like the one 
here considered, the gravitational potential (f) is static 
as long as the density perturbations are linear. Con- 
sequently, in the linear regime the Rees-Sciama ef- 
fect is absent. In our model, the nonlinear regime is 
reached at z ~ 3; nonetheless, the underlying grav- 
itational field is still well described by linear theory 
((/) ^ 1). On sub-horizon scales, this implies that the 
system can still be treated with Newtonian gravity in 
comoving coordinates (Mukhanov et al. 1992). Dur- 
ing the nonlinear epoch, the gravitational field can 
become significantly non-static due to the evolution 
of the density field. CMB photons, which mirror the 
gravitational field (f), can thus refiect not only the sig- 
nature of the gravitational potentials at the surface of 
last scattering, namely the Sachs- Wolfe effect (Sachs 
and Wolfe 1967), but also the imprint from nonlinear 
density perturbations at low redshifts. 

Our technique for calculating CMB temperature 
anisotropies proceeds in two stages (see Anninos et 
al. 1991 for a detailed description of the numerical 
scheme and tests). First, starting at a redshift of z = 
1900, just prior to recombination, a 64^ particle-mesh 
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model is used to simulate the CDM from an initial 
Harrison-Zel'dovic spectrum modified by the appro- 
priate damping term. Since only sub-horizon scales 
are considered, the collisionless matter is evolved us- 
ing Newtonian gravity in an expanding grid. The sim- 
ulation evolves from the linear, through the weakly 
clustered, into the nonlinear regime and terminates 
at z=0. The initial amplitude of the perturbations is 
the only free parameter in the model and was cho- 
sen so that the r.m.s. overdensity in Sh~^Mpc scales, 
(78, reaches unity at z=0. We have verified that this 
normalization yields both reasonable large scale struc- 
ture and peculiar velocity fields for our simulations. 
The results presented here are from CDM models in 
a computational grid 12Ah~^ Mpc on a side, so the 
minimum resolution is ~ 2h~^ Mpc. In the second 
stage of our method, a bundle of photons is propa- 
gated back in time through the evolving dark matter 
structures, starting at at z=0 and finishing at z=1900. 
We propagate bundles with up to 10® individual CMB 
photons, updating their individual temperature at ev- 
ery timestep by computing explicitly the Rees-Sciama 
effect from the second term in equation (1). 

The correlation between the Rees-Sciama effect 
and density inhomogeneities are performed via the 
column density contrast of dark matter along the pho- 
ton trajectories. This column density contrast is de- 
fined as 

. = f '-^dt, (2) 

Jr P 

where the integral is along the photon trajectories. 
Last, we construct maps of AT/T(Rees-Sciama) and 
the column density contrast a. 

Figure 1 shows the r.m.s. value of the Rees-Sciama 
effect as a function of angular scale of the photon 
bundle. The lower limit of the dynamical range of 
the simulation, which is ~ 0.1°, is dictated by the 
finite resolution of our computational model. The 
upper limit of the dynamical range, ~ 8°, is set to 
avoid boundary effects from our periodic computa- 
tional domain. The error bars indicate the cosmic 
variance from 42 simulations, that is, the variance 
in the Rees-Sciama effect for different directions at 
a fixed angular scale. Typically, these variances are 
below the 15% level for photon bundles with angular 
scales > 2°. For a particular photon bundle direction 
on angular scale 2° , figure 2 shows the accumulated 
r.m.s. values of the Rees-Sciama effect {top) and col- 
umn density contrast a {bottom) as a function of red- 



shift. Other directions were considered, and a similar 
outcome was obtained. Three features are evident 
from figure 2. First, as expected, the Rees-Sciama is 
produced almost entirely during the nonlinear stage 
of the density fields (z ^ 3). Second, increases in the 
accumulation of the Rees-Sciama effect are correlated 
with increases in the column density contrast. These 
jumps are due to high density structures within the 
photon beam. Third, because comoving voids make 
the CMB photon appear warmer while high density 
structures make the photon appear colder (Dyer and 
Ip 1988), the accumulation of the Rees-Sciama effect 
is not exactly monotonic; clusters, however, dominate 
the Rees-Sciama signal due their high density con- 
trast {8p/p ^ 10) compared to our ~ 30 Mpc voids 
(0 <^ 8p/ P ^ ~1)) which are still in the quasi-linear 
regime. 

The comparison between the average density con- 
trast 6p/p and the r.m.s. change of the gravitational 
potential ^ within the photon bundle for z < 0.1 
is depicted in figure 3. Here the voids, clusters and 
filaments that the photon beam traverses are directly 
mirrored in the Rees-Sciama signal through We 
see clearly from this figure that indeed overdense 
structures provide the major contribution to the Rees- 
Sciama effect. 

Figures 4 through 7 are maps from the propaga- 
tion of 10® photons as seen by a comoving observer 
at z = 0. Figures 4 and 5 show respectively the Rees- 
Sciama effect and column density contrast a up to 
z = 1900 of a 2° X 2° photon beam. Figure 5 repre- 
sents a map of the total line-of-sight density fiuctu- 
ations, i.e., the galaxy clusters and voids "seen" by 
the CMB photons in figure 4. The poor map corre- 
lation between figures 4 and 5 is due to the smear- 
ing by the void-cluster-void-cluster cumulative Rees- 
Sciama fiuctuations. Figure 6 and 7 show again the 
Rees-Sciama effect and column density contrast a, re- 
spectively. In this case however, the maps are for 
several isolated structures ~ I0h~^ Mpc across us- 
ing a 16° X 16° photon beam. The slab is located at 
0.065 < z < 0.095, about 330/1"! y^^^ ^-^^-^ ^^^-^ ^j^g 
observer and nearly 110h~^ Mpc thick. This is the 
region that yields the maximum density contrast in 
figure 3. We have also examined the region around 
z = 0.03 that produces the maximum and maps 
with similar features were obtained. Since the integra- 
tion to generate figures 6 and 7 is done through a vol- 
ume thin enough, smearing effects are avoided in this 
case. There is a clear correlation between the Rees- 
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Sciama effect (Fig. 6) and the coiumn density con- 
trast (Fig. 7), showing that indeed the Rees-Sciama 
effect is caused by the matter structures. In order to 
quantify the correlation between two given maps, we 
define the correlation strength of maps A and B by 
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(4 + ^5)/2' 



(3) 



where Gij and hij are the corresponding pixel values of 
maps A and B normalized to unity, and N the total 
number of pixels. We found a correlation strength 
^ = .71 between figures 4 and 5. In contrast, figures 
6 and 7 have a stronger correlation of ^ = .93. 

A more careful examination of the maps in figures 
6 and 7 shows that the location of the high density 
structures in figure 7 are not centered at the cold 
(dark) spots in the temperature map (Fig. 6) as one 
would expect. Furthermore, most of the dark spots 
in the temperature map that are located near high 
density structures have associated with them a hot 
(white) spot, with the structure centered between ad- 
jacent hot and cold spots. This configuration resem- 
bles a hot-cold temperature dipole. These dipole-like 
patterns in temperature fiuctuations are due to the 
transverse peculiar bulk motion of the cluster and 
its associated gravitational potential. Photons cross- 
ing the path of the cluster ahead of the approaching 
cluster will be redshifted because of the increasing 
depth of the potential well during the photon cross- 
ing time, while photons crossing at the trailing end 
of the structure are blueshifted. Hence, even an in- 
trinsically static potential can cause a significant tem- 
perature anisotropy as long as the transverse proper 
motion and mass of the cluster are large enough. The 
separation between dipole peaks is indicative of the 
characteristic length-scale of the density perturba- 
tion; furthermore, the temperature difference between 
the dipole peaks is a measure of the transverse veloc- 
ity of the cluster. A detailed analysis of these dipole 
patterns in the CMB sky due to proper motion of the 
inhomogeneities will be presented elsewhere (Laguna 
and Tuluie 1995). 

One can get an estimate of the total effect from 
proper motions and intrinsic potential variations from 
the following scaling argument. The contribution to 
the temperature anisotropy due to an isolated struc- 
ture is 



AT 



5^ 
dt 



6t 



where 6t is the photon crossing time, v is the pe- 
culiar velocity of the structure, d its physical size 
(d ~ 6t) and tc ~ d/(f)^^^ its characteristic evolu- 
tion time (Panek 1992). The first term in equation 
(4) arises from proper motions and the second from 
intrinsic changes in the potential. On the other hand, 
from linear perturbation theory we have that v ~ 
</)(l + z)-3/2(rfF„)-i and 6p/p ~ <j>(l + z)-^(dH,)-\ 
so equation (4) can be rewritten as 



AT 



10" 



3/2' 



rf(i + z)3/2 y 

Uh-^Mpc ) 



+ t)st^vct, + ct,^/', (4) 



(5) 

Within the limitations of this simple argument, we see 
from equation (5) that the effect from proper motions 
is as important as that arising from intrinsic changes 
in the gravitational potential of the structures; more- 
over, if the cluster or supercluster under consideration 
is highly nonlinear, the effect from proper motions 
could have the dominant role in generating these sec- 
ondary perturbations to the primeval CMB. 

In summary, our study considers the gravitational 
infiuence of the evolving density field on the CMB 
directly from the dark matter dynamics and does 
not rely on any assumed symmetry of the struc- 
tures. The underlying matter structures used ex- 
hibit the sheet-filament-knot structure of voids, clus- 
ters and superclusters, providing then a more real- 
istic representation of the shape of density inhomo- 
geneities in the universe. Furthermore, our numeri- 
cal evolution avoids artificial splitting between linear, 
weakly clustered and nonlinear epochs. Our simula- 
tions indicate that the Rees-Sciama effect is small, 
10"^ ^ AT/T ^ 10-^ for void scales ~ 30 - GO/i"! 
Mpc and clusters and superclusters ^ 15/i~^ Mpc. 
Our results disagree with those from authors who use 
larger estimates for the scale and density contrast of 
large scale matter structures (Meszaros 1994; Rees 
and Sciama 1968). For the scales of the structures we 
considered (~ 10h~^ Mpc), we found that the Rees- 
Sciama effect, in addition to being produced by the 
intrinsic evolution of the density field (e.g. changes 
in the depth of the potential wells), is created by the 
change of the potential as the cluster moves across 
the microwave sky. We have found that this effect 
can be as important as that from the intrinsic evo- 
lution of the gravitational potentials in generating 
second-order temperature anisotropies. Finally, our 
simulations show that AT/T anisotropies induced by 
peculiar velocities leave a dipole-like imprint. 
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Fig. 1. — CMB AT/Trms from the Rees-Sciama ef- 
fect as a function of angular scale. The error bars 
are the cosmic variance from 42 simulations using 10® 
photon beams. The minimum and maximum angu- 
lar scales are set by the resolution of the simulations 
~ 2h~^ Mpc and the periodicity of the computational 
domain, respectively. 
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Fig. 2. — Accumulated r.m.s. anisotropy of the Rees- 
Sciama effect (top) and average column density con- 
trast a {bottom) as a function of redshift. The results 
were obtained from a 10® photon bundle with angular 
scale of ~ 2°. As expected, the Rees-Sciama effect in- 
creases significantly during the nonlinear stage of the 
matter evolution z < 3. The contribution from over- 
dense regions is more dominant, so complete cancel- 
lation between dense regions or voids is not observed. 
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Fig. 3. — Comparison of ^ r.m.s. (top) with the 
average density contrast 6p/p {bottom). There is a 
direct correlation between the voids and clusters that 
the photon beam traverses with the rate of change of 
the potentials (-1^). Overdense structures provide the 
major contribution to the Rees-Sciama effect. The 
results were obtained from a 10® photon bundle with 
angular scale of ~ 2°. 



Fig. 4. — Temperature map of 10® photons in a ran- 
domly chosen direction of the simulated microwave 
sky for a 2° x 2° window showing the Rees-Sciama 
effect only. For this and subsequent maps, maximum 
and minimum values are set black (cold) and white 
(hot), respectively. Here, AT/Trms = 2.9 x 10~^ and 
-8.9 X 10"^ < AT/T < 9.5 x 10"^. 
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Fig. 5. — Map of the column density contrast, a, for 
figure 4. This map provides the line-of-sight column 
density contrasts seen by the photons. Here, Uavg = 
-9.5 Mpc and -398 Mpc < (t < 822Mpc. 




Fig. 6. — Map of the Rees-Sciama effect for an iso- 
lated region ~ llO/i-i Mpc thick located at z ~ 0.08, 
containing several typical clusters and voids. The an- 
gular size of the window in this case is 16° x 16°. 
Here, AT/T™, = 7.6 x 10"^ and -2.6 x lO""^ < 
AT/T < 3.7 X 10-"^. Most of the anisotropies are due 
to the dominant clusters seen in figure 7. The tem- 
perature fiuctuations have dipolar features produced 
by the proper motions of the clusters across the sky. 
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Fig. 7. — Column density contrast, a, in dimension- 
less units for the localized region in figure 6. High 
density regions lie at the center of the dipolar pat- 
terns in figure 6. Here, Uavg = 16 Mpc and —70 Mpc 
< (T < 832 Mpc. 



